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ABSTRACT 

The buoyant rise of hot plasma bubbles inflated by AGN outflows in galaxy clusters can heat 
the cluster gas and thereby compensate radiative energy losses of this material. Numerical 
simulations of this effect often show the complete disruption of the bubbles followed by the 
mixing of the bubble material with the surrounding cluster gas due to fluid instabilities on the 
bubble surface. This prediction is inconsistent with the observations of apparently coherent 
bubble structures in clusters. We derive a general description in the linear regime of the growth 
of instabilities on the surface between two fluids under the influence of a gravitational field, 
viscosity, surface tension provided by a magnetic field and relative motion of the two fluids 
with respect to each other. We demonstrate that Kelvin-Helmholtz instabilities are always 
suppressed, if the fluids are viscous. They are also suppressed in the inviscid case for fluids of 
very different mass densities. We show that the effects of shear viscosity as well as a magnetic 
fields in the cluster gas can prevent the growth of Rayleigh-Taylor instabilities on relevant 
scale lengths. R-T instabilities on pc-scales are suppressed even if the kinematic viscosity of 
the cluster gas is reduced by two orders of magnitude compared to the value given by Spitzer 
for a fully ionised, unmagnetised gas. Similarly, magnetic fields exceeding a few fjG result in 
an effective surface tension preventing the disruption of bubbles. For more massive clusters, 
instabilities on the bubble surface grow faster. This may explain the absence of thermal gas 
in the north-west bubble observed in the Perseus cluster compared to the apparently more 
disrupted bubbles in the Virgo cluster. 

Key words: instabilities - galaxies: clusters: individual (Virgo, Perseus) - galaxies: active - 
cooling flows 



1 INTRODUCTION 

Over the last few years numerous models have been proposed to 
explain the observed effects of non-gravitational heating of the 
gas in the centre of galaxy clusters. For a review of the obser- 
vational evidenc e and the a ssociated constraints on heating mod- 
els see Bohringer et al. (2002) and references therein. One of the 
most popular ideas involves the heating of the cluster gas by the 
outflows from AGN at the centres of clusters. The wide variety 
of outflows observed to originate in AGN has led to a number of 
suggestions for the detailed geometry of this model. These range 
from continuous, distributed heating throughout the cluster gas 
through the complete disruption of the outflo w close to the cen tre 
of the cluster I Ruszkowski & Begelmaj.l200llBruggerll2003al) to 
episodic heating by collimated jets comparable to those of powerful 
radio galaxies ( Reyno lds et all l200lt iBasson & AlexandeJ. |2003t 
lOmma et afll2004l) . 

The model geometry in which an episodic AGN inflates bub- 
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bles of relativistic plasma which subsequently detach from the clus- 
ter centre and rise bu oyantly in the c l uster atmosphere, has re- 
ceived most attentionjChAjrazov_et_all l200l t loivflis_etal l_[200lj: 
I Saxton et all l200lb iBruggen & Raised 1200 it brugger! l2003b: 
iDalla Vecchia et all |2004 iRevnolds et all |2004 . There are a 
number of reasons for this concentrated effort. Firstly, although 
the complexity of the problem of modelling outflows from AGN in 
cluste r atmospheres requires numerical simulations (e.g. iBinnevl 
2002), the geometry of individual bubbles rising passively due 
to buoyancy is comparatively simple to investigate. Secondly, al- 
though r adio sources are common at the centres of cooling flow 
clusters <Burnslll99d) . they rarely are powerful enough to inflate 
lobes on scales of several tens to hundreds of kpc. Finally, the ob- 
servation of ellipsoidal depressions of the X-ray surface brightness 
profiles of clusters, sometimes associated with faint, low-frequency 
radio emission suggests this type of inter action of the AGN o utflow 
with the cluster gas (for an overview see lBirzan et allEoo4) . 

A crucial factor in the heating model of buoyant bubbles is 
whether or not the bubbles are stable against break up by fluid insta- 
bilities on their surface. Clearly, in the case of complete disruption 
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of the bubbles and mixing of the bubble material with the surround- 
ing cluster gas, all energy stored in the bubbles can be dissipated. 
While this may be desirable from the point of view of efficiently 
heating the cluster gas, the apparently coherent bubble structures 
observed in clusters argue against complete mixing. Of course, it 
would be difficult in practice to find observational evidence in the 
cluster gas for bubbles after their disruption. We can therefore not 
rule out the possibility that clusters without observational signa- 
ture of buoyant bubbles are still heated by this mechanism. For 
exa mple, the large-sca le radio structure of M87 in the Virgo clus- 
terJOwen^yill [2000]) sugges t s the presence of buoyant bubbles 
I Chur azov et all ' 2001 : Kaiser, 2003), but there appears to be no 
corresponding depressio n in th e X-ray surface brightness profile of 
the cluster (You ng et alll20o3) . 

In this paper we investigate with an analytical approximation 
the stability of the surface of buoyant bubbles rising in cluster at- 
mospheres against disruption by fluid instabilities. We take into 
account the effects of shear viscosity of the cluster gas, which 
was recently suggested as a way to stabilise the bubbles sur- 
face iRevnolds et all 12004V We also study the role of the cluster 
magnetic field by complementing and extending the analysis of 
iDe Young 1 2003). For simplicity in this analytical model, we as- 
sume a static cluster atmosphere in hydrostatic equilibrium in its 
gravitational potential well. There are no large-scale flows or tur- 
bulent motions present in the cluster gas. Such fluid flows may be 
caused by the continuing infall of gas, if the cluster is still in the 
process of forming, or by earlier outflows driven by the AGN. They 
may lead to a faster disruption of the buoyant bubbles, but their 
treatment is beyond the scope of analytical models. 

In Section|2|we discuss the relevant fluid instabilities. Section 
|3|contains a brief discussion of the relevant transport coefficients in 
ionised gases. We apply our estimates to the Virgo cluster in Sec- 
tion |4| and compare our findings with the situation in the Perseus 
cluster in Section [5] We summarise our conclusions in Section [6] 
The appendix contains a detailed derivation of the growth rates of 
fluid instabilities in the linear regime taking into account the ef- 
fects of a gravitational field, surface tension, viscosity and relative 
motion between two fluids. 



2 INSTABILITIES ON THE BUBBLE SURFACE 

Formally the surface of a buoyant bubble in the atmosphere of a 
cluster of galaxies represents a boundary between two fluids. In 
the presence of a gravitational field, this boundary will be subject 
to Rayleigh-Taylor (R-T) instabilities. Since the bubble is moving 
through the cluster atmosphere, Kelvin-Helmholtz (K-H) instabil- 
ities may also develop. Further complications arise from the pos- 
sible presence of a magnetic field in the cluster gas and from the 
potentially relativistic equation of state of the material inside the 
bubbles. 

In the appendix we derive a general expression for the rate of 
growth, n, for instabilities of the boundary between two fluids as a 
function of the wavenumber of the instabilities, k, taking into ac- 
count a gravitational field, viscosity, surface tension and motion of 
the two fluids relative to each other. Here we are mainly interested 
in the situation where the heavier fluid, i.e. the cluster gas, has a sig- 
nificantly higher mass density then the other fluid, i.e. the gas inside 
the bubbles. For this case we find in the appendix that the properties 
of the bubble material are unimportant in determining the growth 
rate of instabilities. Therefore all fluid quantities discussed in the 



following always refer to the heavy fluid, i.e. the cluster gas, unless 
stated otherwise. 

The results in the appendix only apply to the non-relativistic 
case. However, in the simplest case where there is no surface ten- 
sion, no viscosity and no motion of the two fluids relative to each 
other, the growth rate is given by the familiar form n = \/gk, where 
g is the gravitational acceleration at the fluid boundar y. In the rela - 
tivistic case this expression is modified to I Allen & Hughes! 1 1984 



gkp 

Sp/c 2 + p' 



(1) 



where p and p are the rest mass density and the pressure of the 
cluster gas, respectively. For typical values of the particle number 
densities of 10~ 2 cm~ 3 and pressures of 10~" ergscm~ 3 in the 
cluster gas, the relativistic correction term 8p/c 2 , equivalent to the 
relativistic inertia of pressure, is at least five orders of magnitude 
smaller than p. In our regime we can thus safely neglect any correc- 
tions arising from a relativistic treatment of the gas in the bubbles. 

In the appendix we also show that for fluids of very differ- 
ent mass densities K-H instabilities are suppressed completely on 
the boundary between the fluids. This result also holds for inviscid 
fluids. Therefore we will neglect K-H instabilities in the analysis 
below. 

For practical purposes the wavenumber, k, can be transformed 
into the corresponding length scale of the perturbation, / = 2n/k. 
The growth rate n gives the growth time f grow = \/n which mea- 
sures the time it takes the amplitude of a given perturbation to grow 
by a factor e. To first approximation the perturbation changes to the 
non-linear regime at this point and grows much faster afterwards. 
Howe ver, the inst ability takes a few e-folding times to fully develop 
JPe YounglEool) . 



3 TRANSPORT COEFFICIENTS IN THE CLUSTER 
ATMOSPHERE 

The temperature in clusters is a few to several keV, implying fully 
ionised g as. The visco sity of ionised, unmagnetised hydrogen is 
given by <SpitzeitlT962T) 



,u = 0.406 - 



e 4 lnA scm 



(2) 



where ra p is the proton mass, k^ is the Boltzmann constant, T is the 
temperature of the gas and e is the elementary charge. For our case 
of a pure hydrogen gas the Coulomb logarithm is 



A = 247trc e 
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where n e is the electron number density of the gas. The kinematic 
viscosity appearing in the equations of the appendix is then v = 

v/p- 

The importance of viscous processes in the transport of mo- 
mentum in a fluid compared to inertial effects on typical length 
scales L is given by the dimensionless Reynolds number, Re = 
uL/v. u is the typical velocity of the fluid on scales comparable 
to L. For cluster gas of a temperature k^ T ~ 3 keV and an electron 
density n e ~ 0.1 cm~ 3 , we find Re ~ 100 for transonic velocities 
on scales of roughly 1 kpc. Obviously, Re will be even smaller for 
subsonic motion. We therefore conclude that viscosity can play an 
important role in the fl uid dynamics of the cluster gas ( Fabi an et all 
120031 iRuszkowski et allE004l) . if the kinematic viscosity is close 
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to the value derived by Spitzer, i.e. if the cluster gas is fully ionised, 
but does not contain a magnetic field. 

Thermal conduction is not directly relevant for the stability or 
otherwise of the surface of buoyant bubbles in clusters discussed in 
this paper. Nevertheless, it is interesting to note here that the trans- 
port of heat by conduction in the cluster gas is intimately linked 
with viscosity as both depend essentially on the same microscopic 
processes between gas particles. In t he case of an unmagnetised 
plasma as described by Spitzer 1 1962), the thermal conductivity of 
the plasma is given by 



= 4.5 



3/2 



(k B T) 5 ^ 2 k B ergs 



/mTe 4 lnA scmK 



(4) 



where m e is the mass of an electron. The thermal diffusivity of the 
plasma is % = k/ (t' p p) , with 



5k B 
2// m m p ' 



(5) 



the specific heat at constant pressure. /j m is the mean mass of a 
gas particle in units of m p and for pure ionised hydrogen /j m ~ 
0.6. To compare the relative importance of viscosity and thermal 
conduction for the transport of energy in the gas we can form the 
dimensionless Prandtl number Pr = v/%. Using equations and 
we find 



Pr~0.4, 



'0.01. 



(6) 



To first approximation Pr does not depend on the exact hydrody- 
namical conditions within the plasma, at least if the plasma is un- 
magnetised. This result means that if viscous effects are important 
in the cluster gas, then thermal conduction of heat will also play a 
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the centres of cooling flow clusters. If, for example, the viscous dis- 
sipation of the energy stored in sound waves excited by the outflows 
from AGN contributes t o the heating of cluster gas l Fabia n et all 
l2003HRuszkowski et all E004). then thermal conduction should be 
important too. However, because Pr is small, this line of reason- 
ing cannot be reversed. The gas in cluster atmospheres, in which 
thermal conduction is important, is not necessarily viscous. 

In the case of a magnetic field in the plasma, the thermal con- 
duction of heat and therefore also the viscous transport of momen- 
tum is suppressed perpendicular to the fieldlines. In a magnetic field 
tangled on scales much smaller than the typical size of the plasma 
th is leads to a reduced value for t he average ther mal condu ctivity, 
V JChandran & Cowlevlll99l:lNaravan & Medvedevll200ll) . How- 
ever, in the presence of turbulent motions in the clu ster gas, anoma- 
lous heat transport may also occur ICho et aill2003l) . increasing the 
effective thermal conductivity of the gas. Similar effects are likely 
to affect the viscosity of the gas. It is not clear which description 
of thermal conductivity or viscous transport of momentum is ap- 
plicable to the magnetised plas ma in galaxy c l usters (for the issue 
of thermal conductivity see Voiat et all l2002t iFabian et all 12002b 
IVoiet & Fabiari |2004) and the unique conditions in individual clus- 
ters may well rule out a generic answer to this question. 



4 APPLICATION TO THE VIRGO CLUSTER 

We now proceed to apply the results of the appendix and Section 
[3] to the hydrodynamical conditions derived from X-ray observa- 
tions in galaxy clusters. In most clusters the buoyant bubbles are 
observed as depressions in the X-ray surface brightness of the clus- 
ter gas and are locat ed close t o the centres, ~ 20kpc, of the re- 
spective clusters iBirzan et alll2004l) . Therefore we need reliable 
estimates of the gas density and temperature in the cluster atmo- 
sphere down to very small radii. This requirement is aggravated by 
the need to determine the gravitational acceleration of the cluster 
gas from the spatial derivatives of these properties. 



4.1 Observational constraints 

The high-resolution X-ray maps of the nearby Virgo cluster pro- 
duced by CHANDRA and XMM-NEWTON allow the determination 
of the gas properties down to a few kpc from the cluster centre 
(Matsushita et all l2002tfYoung et allEoollGhizzardi et alll2004 
hereafter G04). G04 obtain a good fit to the electron density distri- 
bution using the fitting function 
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withw! = 0.089 cm~ 3 ,f; 2 = 0.019 cm" 3 , r x =4.5 kpc, r 2 = 21 kpc, 
a; = 1.5 and 0C2 = 0.7. For the temperature distribution they use 

2 J. 



r = T - 7i exp I ---2 



with k B T = 2.4keV and k B T { = 0.78 keV. 

Assuming hydrostatic equilibrium for the gas in the gravita- 
tional potential well defined by the dark matter halo of the Virgo 
cluster, it is also possible to determine the mass distribution (G04) 
and therefore the gravitational acceleration g within this cluster. 
The result from this analysis agrees very well in the region of 
overlap with the mass estimates derived from the dy namics of the 
globu lar clusters and the stars of M87 iRo manowskv & Koc hanek. 
2001). The well known properties of the Virgo cluster make it an 
ideal example for the analysis of the stabili ty of a buoyant bubble. 

The torus-shaped radio structure lOwen et all l2(X)ro associ- 
ated with the central galaxy M87 on scales of approximately 15 kpc 
to the east of the cluster centre has been interpreted as a buoyant 
bubble caused by the AGN at its centre IChurazov et all 1200 it 
iKaiseil 120031) . The situation is less clear to the west of M87, 
but a similar, more disrupted buoyant structure may be present 
there. However, no corresponding depressions of the X-ray surface 
brightness have been observed lYoung et al., 2002). The eastern ra- 
dio structure resembles a torus, inclined towards our line of sight. 
This is fully consistent with the structure expected for a buoyant 
bubble deformed, but not yet fully disrupted, by fluid instabilities 
on scales comparable to the bubble size of ~ 10 kpc IZhidov et all 
1977; Churazov et al., 2001). On smaller scales, below the resolu- 
tion limit of current observations, fluid instabilities may already be 
fully developed and may have mixed the bubble material with the 
cluster gas. The relativistic, magnetised plasma initially injected 
into the bubbles gives rise to the synchrotron emission while the 
mixed-in thermal cluster gas is responsible for the X-ray emis- 
sion observed in the same region. The elongated radio emission 
region connecting the buoyant bubbles with the centre of M87 may 
be caused by small pockets of relativistic plasma created and de- 
tached from the main bubbles by fluid instabilities operating on 
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small scales. We would not expect to observe these features 'trail- 
ing' behind the main bubbles if the bubbles were still fully intact. 
Thus the buoyant bubbles in Virgo represent a good test case for 
our study as they suggest that fluid instabilities on the surface of 
the currently rising buoyant bubbles have fully developed on pc- 
scales, but are only starting to disrupt the bubbles on kpc-scales. 

The disruption and mixing of the buoyant bubbles with the 
cluster gas may influence the observed radio emission. The bub- 
bles were originally inflated by the central AGN with a magnetised 
relativistic plasma similar to that of the currently 'active' lobes on 
scales of ~ 3 kpc surrounding the observed jet and postulated coun- 
terjet. We would expect the currently observed synchrotron radio 
emission to originate from the initial bubble material and not from 
any mixed-in thermal cluster gas. If the mixed-in cluster gas was 
also strongly magnetised, then internal depolarisation should lead 
to very low level s of polarisa tion of the radio emission from the 
buoyant bubbles l Burg, Il96a) . The fact that the synchrotron ra- 
dio emission from the radio structures interpreted here as buoy- 
ant bubbles is highly polarised at high frequency (10.55 GHz, 
iRottmann et all [l996h . seems to argue against the bubbles being 
well mixed in with the cluster gas. However, the interpretation of 
synchrotron radio polarisation data depend s strongly on the exact 
assumed geometry of the emitting region lLainsl ll984l) . Further- 
more, if the mixed-in cluster gas is not or only weakly magnetised, 
then it will not have a significant effect on the polarisation of the 
synchrotron emission from the bubble material. We will see below 
that our findings from the stability analysis of the bubble surface ar- 
gue for at best a weak magnetic field in the cluster gas surrounding 
the bubbles. 



4.2 The effects of viscosity and a magnetic field 

The solutions of equation <A32> clearly depend on the location 
within the cluster under consideration as all relevant cluster prop- 
erties change as a function of distance r from the cluster centre. 
This implies that instabilities on the bubble surface of different 
wavenumber k, or, equivalently, different length scale I = 2%/k, 
will be able to grow at different r. The buoyant bubbles are rising in 
the cluster atmosphere and so instabilities may start growing when 
a given bubble is close to the cluster centre, but may be 'frozen' 
again as the bubble travels further out. Alternatively, instabilities of 
a certain scale length may be suppressed in the cluster centre, but 
start growing once the bubble rises to larger r. 

In the following we use the analytical fits to the electron den- 
sity and gas temperature in the Virgo cluster presented by G04, 
equations Q and {SJ. Figure Q shows the resulting gravitational 
acceleration and the kinematic viscosity for an unmagnetised hy- 
drogen plasma in hydrostatic equilibrium within the Virgo cluster. 
Using these results we now estimate the growth time, ? gl - ow , for fluid 
instabilities from equations JA32> . <A35> and <A36> as a function 
of their size /. Figure [2] illustrates the effects of the viscosity and 
a magnetic field on the growth time of the instabilities of various 
sizes. Note that these solutions only apply to a bubble located at 
a distance of 10 kpc from the centre of the Virgo cluster. The sit- 
uation changes considerably with the bubble position within the 
cluster (see below). 

In the absence of viscosity and surface tension from a mag- 
netic field, instabilities grow on all size scales with the familiar 
proportionality of f grow \fl. Viscosity and magnetic fields only 
affect the growth of instabilities on small scales. The introduction 
of a magnetic field in the cl uster gas results i n a sharp cut-off at 
small sizes (see also C61 and lDe Youngll2003l) . For illustration we 
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Figure 1. Gravitational acceleration (solid line) and kinematic viscosity 
(dashed line) of the gas in the Virgo cluster. Both quantities are derived 
from the analytical fit to the electron density and temperature profiles based 
on deprojected X-ray CHANDRA data by G04, equations an£ l 111- The 
vertical dotted lines show the position of the data points used by G04. The 
kinematic viscosity is calculated for the case of a fully ionised, unmagne- 
tised hydrogen plasma, equation 01 . 
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Figure 2. Growth time of fluid instabilities on the surface of a buoyant 
bubble at a distance of 10 kpc from the centre of the Virgo cluster. Solid 
line: Real part of the fastest growing solution of equation IA32I taking into 
account the viscosity of the cluster gas and a magnetic field of strength 1 fjG. 
Dotted line: Same as solid line, but neglecting the surface tension provided 
by the magnetic field. Short-dashed line: Same as solid line, but neglecting 
viscosity. This is the solution of equation IA35I . Long-dashed line: Growth 
time of instabilities without viscosity or magnetic field. Solution to equation 
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used a magnetic field of 1 ,uG in Figure|2| Here and in the following 
we assume that the magnetic field is parallel to the bubble surface. 
This assumption is justified as any field components normal to the 
bubble surface will not influence the growth of instabilities in the 
linear regime (C61). Note that even such a weak magnetic field has 
a profound effect on the growth of the instabilities. Our assumption 
that the density of the cluster gas strongly exceeds that of the gas 
in the bubble leads to the disappearance of all quantities associated 
with the bubble material from the relevant expressions (see the ap- 
pendix). The same holds for the magnetic field. Thus the magnetic 
field discussed here is the field in the cluster gas, not the field inside 
the bubble. The latter may be by far the stronger of the two, but, in 
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our approximation, does not influence the growth of instabilities on 
the bubble surface. 

The cut-off due to the magnetic surface tension below which 
instabilities cannot grow is given by setting fgiow — * 00 or, equiva- 
lently, n = in equation <A35t . Thus 



and we would not expect a given bubble to disintegrate into units 
smaller than typically / m ; n (see also lDe Yound,l2003l) . The break up 
of the bubble into smaller bubbles of sizes only slightly larger than 
Z m j n happens rapidly. In fact, the minimum growth time is reached 
for an / only slightly larger than / m ; n . The structure of the bubble 
on large scales remains coherent for longer. If the bubble were ob- 
served from a distance with limited resolution, it would give the im- 
pression of being intact, while on small scales the observed large- 
scale structure is made up of many smaller units with a matrix of 
intermixed cluster gas in between them. Similar structures result 
when a coh erent bubble of air is suddenly released under water 
IZhidov et ail 1 19771) . 

In the case of a viscous gas without a magnetic field, the 
growth of instabilities on the bubble surface is significantly delayed 
on small scales. There is no cut-off as in the case of a magnetic 
field, but again the growth time reaches a minimum and no insta- 
bilities on any size scale can grow faster than this minimum time. 
Finally, in the case of a viscous gas containing a magnetic field, 
the sharp cut-off introduced by the surface tension provided by the 
magnetic field is smoothed out somewhat by the effects of viscos- 
ity. 
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Figure 3. Limits for the scale size, I, of instabilities able to develop on the 
surface of a buoyant bubble in the Virgo cluster as a function of radius 
within the cluster atmosphere. Long-dashed line: Upper limit, l up , from 
equation <A36I neglecting viscosity and magnetic fields. This upper limit 
also applies to good approximation to all other cases considered here (see 
text). Short-dashed line: Lower limit, /i ow , for a viscous, unmagnetised clus- 
ter gas with kinematic viscosities equal to the Spitzer value (top) and re- 
duced by a factor 0.01 (bottom). Solid lines: Lower limit /[ ow for the case 
of an inviscid, but magnetised cluster with field strengths 0.1, 1 and 10 /jG 
from bottom to top. Instabilities on size scales between the lines giving 
lower limits, i[ ow , and the long-dashed line at any fixed distance from the 
cluster centre, r, are able to develop within the rise time of the buoyant 
bubble currently located at r. 



4.3 Stability of buoyant bubbles 

Figure[2]compares the effects of viscosity and the magnetic field in 
the cluster gas on the stability of the surface of a buoyant bubble. 
However, it only provides us with a snapshot at a single time during 
the rise of the bubble. As the bubble rises, the conditions for the 
growth of instabilities change and so we cannot use Figure [2] to 
estimate the degree of disruption of the bubble at a given time. 

In the following we make the simplifying assumption that the 
buoyant bubbl es in the Virgo cluste r rise at a constant velocity 
of WOkms-'. IChurazov et alJfeOOlh showed that their numerical 
simulations of the bubbles in the Virgo cluster are consistent with 
a rise velocity of order the Keplerian velocity, Vr = \fgr, where r 
is the distance from the cluster centre. For the inner 20 kpc of the 
Virgo cluster, is almost constant at 400 km s~ 1 . The position of a 
buoyant bubble starting its rise at t = and r = r\ = \ kpc is simply 
given by r = v%t + r x . We obtain a rough estimate for the wavenum- 
bers of instabilities that are able to grow on the bubble surface by 
setting ? glow = t and solving equations <A32t . <A35> and <A36> for 
k. We then convert the wavenumber to the relevant size scale, /, of 
the instabilities. We will assume for simplicity that the strength of 
the magnetic field in the cluster gas does not change as a func tion 
of distance from the cluster centre (but see lMurgia et alll2004) . 

In the case of an inviscid cluster gas without a magnetic field, 
equation <A36> . this prescription results in an upper limit, Z up , for 
the size of instabilities that are able to develop. Instabilities with 
I < Z up should be well developed in this case. Taking into account 
viscosity or a magnetic field in the cluster gas results in a minimum 
for the growth time as a function of wavenumber or instability scale 
size (see Figure|2j. Depending on the rise time of the bubble, insta- 
bilities between a lower limit /] ow and an upper limit / up can develop 



or the bubble is stable on all scales. Unless the rise time of the bub- 
ble is close to the minimum growth time, Z up is almost identical for 
all cases. 

Figure|3|shows the results for this rough estimate. In all cases 
the upper limit / llp starts to exceed the observed bubble size of about 
10 kpc soon after the bubble starts rising. We therefore expect that 
the bubble should be deformed on scales comparable to its overall 
size. The observed, well-developed torus shape of the eastern bub- 
ble and the apparently more advanced disruption on large scales 
of the western bubble are consistent with this prediction. For small 
instabilities the situation is more complicated. We pointed out in 
Sec tion l4~T1 that the absence of X-ray depressions in the region of 
the radio structures suggests the buoyant bubbles to be well mixed 
with the thermal cluster gas on pc-scales. If the cluster gas has a 
viscosity as high as the value given by Spitzer' s treatment, then no 
instabilities smaller than about lOOpc can develop on the bubble 
surface at any stage during the rise of the bubble. Further out in 
the cluster atmosphere, viscous effects will suppress instabilities 
even up to 1 kpc. For instabilities to be able to grow on pc-scales 
the kinematic viscosity has to be reduced by about two orders of 
magnitude. While such a significant suppression can be expected 
in the case of the related process of thermal conductivity due to the 
prese nce of magnetic fields in the cluster gas IChandran & Cowlevl 
1998), it is currently not clear whether viscous effects are similarly 
reduced. 

Turning to the effects of the magnetic surface tension, we find 
that even weak magnetic fields of IjjG already stabilise the bubble 
surface on scales down to about 10 pc against R-T instabilities. The 
strength of the magnetic field in the cluster gas must be as low as 
about 0. l,uG for instabilities on sub-pc scales to grow. Only for such 
low field strengths can we expect significant mixing of the bubble 
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material with the thermal cluster gas. We will discuss our findings 
on viscosity and magnetic fields in more detail in Section|6| 



5 COMPARISON WITH THE PERSEUS CLUSTER 
5.1 Observational constraints 

The Perseus cluster is significantly further away from us than the 
Virgo cluster. Thus it is impossible to achieve the same accuracy in 
determining the gas properties and gravitational acceleration for the 
cluster regions inside 20 kpc. However, we are mainly interested in 
studying the changes to the stability analysis for buoyant bubbles 
caused by a shallower gravitational potential. Furthermore, while 
we argued that the buoyant bubbles in the Virgo cluster are well 
mixed with the thermal cluster gas, the situation is quite different 
in the Perseus cluster. 

Two bubbles associated with the currently active AGN in the 
central galaxy NGC 1275 are detected both as depressions in the 
X-ray surface brig htness and in synchrotron radio emission (e.g. 
iFabian et alll200(j) . There is at least one other depression in the X- 
ray surface brightness about 15 kpc to the north-west of the cluster 
centre connected by a s pur of low-fre quency radio emission with 
the central radio source (Fabian et al., 2002). In the following we 
concentrate on this outer bubble, since it is located at a similar dis- 
tance from the cluster centre as the buoyant bubbles in the Virgo 
cluster. The 'hole' apparent in the X-ray maps is consistent with 
the responsible buoyant bubble being devoid of any thermal gas 
(Schmidt et al., 2002; Sanders et al., 2004). It is therefore unlikely 
that small, pc-scale instabilities have mixed the bubble material 
with the cluster gas like for the bubbles in the Virgo cluster. On 
larger scales the bubble appears flattened in the rise direction and 
may also have been distorted into a torus shape similar to the east- 
ern bubble in the Virgo cluster. This morphology indicates that the 
north-western bubble in the Perseus cluster may be in the process 
of deforming on scales comparable to the bubble size, ~ 10 kpc, 
while its surface is stable on smaller scales. 

In order to derive the properties of the cluster gas and the grav- 
itational acceleration of the Perseus cluster we use the mean depro- 
jected elec tron dens i ty and temperature profiles presented in Fig. 
21 of ISanders et al] {2004). We fit the data points with the same 
expressions for the electron density and temperature we used for 
the Virgo cluster, equations Q an d {8}. However, there are only 3 
data points inside a radius of 20 kpc and so we simplify the density 
profile by setting »j = 0. We find a good fit for n\ = 0.07 cm~ 3 , 
r\ = 30 kpc and 04 =0.85. The corresponding fit for the tempera- 
ture profile is k% Tq = 14.7 keV, k%T\ = 11.6 keV and ri = 1 14 kpc. 
Note that the temperature of the gas in the Perseus cluster is still 
rising outwards at distances beyond 120 kpc indicating a shallower 
gravitational potential than that of the Virgo cluster. The fits are 
based on only three data points inside the position of the buoyant 
bubble at about 20 kpc from the centre. Clearly this rules out any 
attempt to derive results as accurate as for the Virgo cluster. How- 
ever, we can still use these fits to illustrate systematic changes of 
our analysis for the given different properties of the two clusters. 

From the density and temperature distributions in the Perseus 
cluster we derive the kinematic viscosity and gravitational acceler- 
ation, see Figure|4| The values for both quantities can be compared 
with their values for the Virgo cluster, FigureQ The kinematic vis- 
cosity is comparable in both clusters, but the gravitational acceler- 
ation is much lower at a given radius in the Perseus cluster. 
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Figure 4. Gravitational acceleration (solid line) and kinematic viscosity 
(dashed line) of the gas in the Perseus cluster analogous to Figure HI for 
the Virgo cluster. The vertical dotted lines show the position of the data 
points inside 20 kpc presented by Sanders e t all l2004h . 
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Figure 5. Limits for the scale size, I, of instabilities able to develop on the 
surface of a buoyant bubble in the Perseus cluster as a function of radius 
within the cluster atmosphere. Lines are as in Figure |3] but calculated for 
the conditions in the Perseus cluster. 

5.2 Stability of buoyant bubbles 

Due to the monotonic behaviour of the gravitational potential inside 
20 kpc, the Keplerian velocity is monotonically rising. However, 
given the uncertainty for the value of the gravitational acceleration 
and for simplicity we make the same assumption as for the Virgo 
cluster that the rise velocity of any buoyant bubble is constant. We 
set it to 300 km s -1 . 

FigurefJJshows the limits on the sizes of fluid instabilities that 
can grow on the surface of a buoyant bubble in the Perseus clus- 
ter. This is directly comparable to the situation in the Virgo clus- 
ter presented in Figure [3] The upper size limit for the instabilities 
able to grow is lower in the Perseus cluster compared to the Virgo 
cluster. However, we would still expect well-developed instabilities 
on scales of several kpc at radii ~ 20 kpc, consistent with the ob- 
servations of the north-western bubble in the Perseus cluster. Vis- 
cous effects prevent the growth of instabilities below a size of about 
1 kpc if the kinematic viscosity is equal to that of an unmagnetised 
plasma. Even for a viscosity reduced by two orders of magnitude 
the bubble surface is stable down to several 10s of pc. The apparent 
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absence of thermal material from the bubble on small spatial scales 
could therefore be explained by viscous effects. However, the sur- 
face tension provided by the cluster magnetic field would have the 
same effect for field strengths as low as a few tenths of /jG. Mix- 
ing of the bubble material with the cluster gas on sub-pc scales is 
therefore unlikely within the Perseus cluster. This result contrasts 
with the situation in the Virgo cluster and is mainly caused by the 
smaller gravitational acceleration in Perseus. 



6 DISCUSSION AND CONCLUSIONS 

We used the formalism of C61 to derive the linear growth rates of 
R-T and K-H instabilities on the boundary between two fluids. We 
demonstrated that the K-H instability is always suppressed if the 
fluids are viscous. Furthermore, the K-H instability is also negligi- 
ble on all length scales for inviscid fluids, if the heavier fluid has 
a significantly larger mass density than the light fluid. We use our 
linear analysis to construct a simple prescription for the growth of 
fluid instabilities on the surface of light bubbles buoyantly rising 
in galaxy cluster atmospheres. Application of our prescription to 
the conditions derived from X-ray observations in the Virgo and 
Perseus clusters leads us to the following main conclusions. 

In the absence of magnetic fields and highly suppressed shear 
viscosity the surface of buoyant bubbles is highly susceptible to dis- 
ruption by the R-T instability on all size scales up to the size of the 
bubbles. The complete mixing of bubbles with the thermal cluster 
gas will occur before the bubbles can rise beyond the innermost re- 
gions of the cluster. This behaviour is confirmed by numerical sim- 
ulations down to the numerical resolution (e.g. Briieaen & Kaiser, 
2002). The mixing is complete on all length scales such that no 
trace of the initial geomtery of the bubble remains. 

Both shear viscosity and the presence of a magnetic field in 
the cluster gas can prevent the disruption of the bubbles on small 
scales. Viscous effects delay the growth of the instabilities while 
the magnetic field provides an effective surface tension. Mixing is 
prevented even if the kinematic viscosity is significantly suppressed 
by two orders of magnitude compared to the value given by Spitzer 
for a fully ionised gas. In the case of a magnetised cluster gas, a 
field strength of a few /jG is sufficient to stabilise the bubble sur- 
face. 

An important factor in determining whether buoyant bubbles 
are disrupted in a given cluster atmosphere is the strength of the 
gravitational field. The gravitational acceleration sets the lower 
limit for fluid instabilities which are able to grow in the presence of 
a magnetic field via equation The effect of shear viscosity also 
depends on the gravitational acceleration. In general smaller insta- 
bilities can grow faster for larger gravitational acceleration. Thus 
the complete mixing of buoyant bubbles with the cluster gas should 
occur preferentially in cluster with a deeper gravitational potential 
well. 

In the case of the buoyant bubble observed in the Perseus clus- 
ter to the north-west of the cluster centre, the above considerations 
explain the coherent appearance of and the apparent absence of any 
thermal gas inside the bubble ISchmidt et aUl2002tlSanders et all 
2004). The buoyant structures in the Virgo cluster observed in the 
radio, however, are not associated with depressions in the X-ray 
surface brightness and may therefore contain substantial amounts 
of thermal cluster gas. This observation suggests that the shear vis- 
cosity in the Virgo cluster must be highly suppressed to at least 1% 
of the Spitzer value. Furthermore, the magnetic field strength in the 
cluster gas can be at most of order a few 0.1 /jG. 



We cannot measure the kinematic viscosity of the cluster gas 
directly and theoretical efforts to calculate transport coefficients in 
ionised, magnetised gases have so far concentrated on thermal con- 
ductivity. It is thus hard to assess whether a suppression of the kine- 
matic viscosity in the Virgo cluster by at least two orders of magni- 
tude is realistic. Even in the case of thermal conductivity estimates 
for th e suppr ession factor range from 0.00 1 IChandran & Cowlevl 
1998) to 0.4 jNaravan & Medvedev[ [2001). The recent numerical 
simulations by ( Reynolds et al., 2004) suggest that shear viscosity 
in the cluster gas leads to bubble morphologies and flow patterns 
consistent with observations. However, more effort is needed in the 
theoretical treatment of shear viscosity in magnetised plasmas to 
decide its importance in the cluster gas. 

We can estimate the strength of the magnetic field in the 
cluster from rotation measurements, RM, towards the central ra- 
dio source in M87. lOwen et"al] ll990l) find an average RM of 
about 1000 radm~ 2 with peaks around 8000 radm~ 2 . If the Fara- 
day screen responsible is distributed throughout the cluster gas, the 
implied strength of the magnetic field is 40 /jG. However, our line of 
sight to the inner radio lobes of M87 passes through the large-scale 
radio structure which will also act as a Faraday screen. It is possi- 
ble that the entire observed RM originates in the outer radio emit- 
ting regions IChurazov et allEoOll) . In this case, the magnetic field 
in the cluster gas itself is substantially lower. In general various 
methods employed in determining the strength of mag netic fields in 
clusters give different results l Carilli ^Javtoi|2002) and they may 
vary within individual clusters (Murgia et al., 2004). However, it is 
unlikely that the strength of the magnetic field in the vicinity of the 
buoyant bubbles in the Virgo cluster is much less than 1 ,uG. This 
value may just be low enough to allow the mixing of the buoyant 
bubbles with the cluster gas on small length scales. 

In summary we find that buoyant bubbles should in general 
remain intact in cluster atmospheres. Exceptions like the bubbles 
in the Virgo cluster require a suppresion of the shear viscosity by 
more than two orders of magnitude compared to the value given by 
Spitzer. Furthermore, in such cases the strength of the cluster mag- 
netic field must be below a few /jG. It is therefore unlikely that we 
are failing to detect the remains of bubbles in many clusters of mod- 
erate mass because they have been disrupted and mixed in with the 
cluster gas. In the most massive systems with weak magnetic fields 
the disruption of bubbles may occur more easily. In these clusters 
the absence of depressions in their X-ray surface brightness profiles 
does not preclude the possibility that their atmospheres have been 
heated by buoyant bubbles in the past. 

This possibility is significantly enhanced if large-scale gas 
flows or turbulent motions are present in the cluster atmosphere. 
Both types of flows may be set up by gas infall in the earlier stages 
of cluster formation as well as the passage of bubbles caused by 
earlier activity cycles of the AGN through the cluster gas. As an- 
alytical models are unable to account for these additional effects, 
future numerical simulations of buoyant bubbles in clusters must 
address these issues. 
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APPENDIX A: INSTABILITIES AT FLUID BOUNDARIES 



Here we present a combined analysis taking into account the R-T and K-H instabilities as well as the effect s of viscos i ty and surface tension 
at the boundary between two fluids superposed in a gravitational field. We hereby extend the work offchandrase khaJ ll96ll hereafter C61), 
who neglected viscosity in his treatment of the K-H instabi lity. A similar exte nsion also, including a possible rotation of the two fluids, was 
attempted previously by Bhatia & Sankhla 1 1984) and Mehta & Bhatia 1 1988). However, we show below that their treatment of the problem 
is flawed. 

Consider a fluid of mass density p and pressure p in equilibrium in a gravitational field acting along the z-axis of the chosen coordinate 
system. The fluid has a coefficient of shear viscosity of p. For simplicity we assume that p is a function of z only. We assume that the velocity 
field in the fluid remains divergence free at all times, V ■ v = 0, which is equivalent to the fluid being incompressible. This condition is a good 
approximation as long as all motions remain subsonic (e.g. Ishvlll99ll) . The fluid may stream uniformly in the direction of the x-axis at a 
velocity U. To allow for a situation where two fluids are superposed in the gravitational field, all quantities describing the fluid may change 
discontinuously at some xy-plane with coordinate z s - On this surface a possible surface tension T s must be taken into account as well. 

Now consider a small disturbance of the fluid such that its density in the perturbed state is given by p + 8p. The associated pressure 
perturbation is given by hp and the perturbed velocities in the fluid are U + u in the x-direction, v in the y-direction and w in the z-direction. 
We assume that the perturbation is small so that u, v, w <C U. The surface at which discontinuities in the fluid properties occur is displaced in 
the z-direction by a 8z s due to the perturbation. Note here that the velocities u, v and w are perturbations and therefore are continuous at z s - 
The equations governing the subsequent evolution of the perturbation are then given by (check for equations 8 to 12 and 28 in chapter X and 
equations 1 to 6 in chapter XI of C61) 
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where quantities with a subscript V are evaluated at the surface located at Zs- 

We now consider perturbations of all fluid quantities of the form /(z)exp (ik x x + ik y y + nt) , where /(z) is an undetermined function 
which may be different for different quantities. Clearly only perturbations for which n has a positive real part will grow in time. In the 
following we will derive expressions for n. 

Substituting this form for the perturbation into the above equations yields 
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where k = ^Jk x +k 2 ,. By multiplying equation <A7t by — ik x , equation <A8> by —iky, summing and making use of equation <A12> we arrive 
at 



,2c / dv " ,l\ dw dp f d 2 ,,\ ., dU 

k 2 8 p=-n'p- + p(^-k 2 \- + -(^ + k 2 )w + l k x pw-, 



(A13) 



where n' = n + ik x U. Substituting into equation <A9> from equation <A10> for 8p and from equatuon <A1 H for 8z s yields 



d s _l ( d2 ,A .-.dpdw gwdp j2„w s 

-8p = -npw + p(^-k w + 2-- + ---kT s - F 8(z-z s ) 



(A14) 
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We can now eliminate 8p between equations <A13> and <A14> to give 
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Equation IA15I is a differential equation describing the possible solutions for the velocity of the perturbation in the z-direction, w. 
We are not interested in the exact form of w, but the solutions of equation <A15t combined with its boundary conditions allow us to find 
expressions for n. 

Here we are interested in cases where within a gravitational field a heavy fluid with mass density P2 and viscosity fj 2 is placed on top of 
a lighter fluid with mass density pi and /ji . This situation can give rise to the R-T instability. We choose the location of the boundary between 
the two fluids such that z s = 0, so that the heavy fluid is located at z > and the light fluid at z < 0. For simplicity we assume that all fluid 
properties are uniform within the domain of each fluid and only change discontinuously at the boundary. This assumption also implies that 
both fluids move along the boundary at velocities U\ and JJi, respectively. Since in general U\ / U 2 , the resulting shear at the boundary may 
cause K-H instabilities to grow. Note here that n' also changes its value at the boundary. With these assumptions equation IA15i reduces to 
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inside the fluids, away from the boundary. By setting 

q 2 =k 2 

where v = /j/p is the kinematic viscosity, the differential equation is simplified to 
d2 -*V£-/U = 0. 
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The solutions to this equation are given by linear combinations of the terms exp(±fe) and exp(±gz). To avoid unphysical solutions of the 
form w — > <x> for z — > the solutions must have the form 
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Here and in the following subscripts denote quantities appropriate for the two individual fluids. The factors n\ and n' 2 are included in the 
solutions for convenience when applying the boundary conditions at the surface between the fluids further on. 

C61 shows that four boundary condi tions must hold a t the surface between the tw o fluids. Equation <A1U implies that w/n' must be 
continuous at the boundary. Both Bhatia & Sankhla ( 1984) and Mehta & Bhatia 1 1988) apply an incorrect boundary condition (w\ = w 2 ) 
at this point. The same continuity at the boundary applies also to dw/dz because of equation <A12> combined with the continuity of u and 
v. Note that if dw/dz is not continuous or ill-defined at the boundary, then the last term of equation <A15> implies that viscous stresses at 
this boundary grow without bounds until the continuity of u, v and dw/dz is established. C61 shows that the continuity of viscous stresses 
tangential to the boundary also imply that fj (d 2 /dz 2 + k 2 } w is continuous across the boundary. The final condition arises from integrating 
equation JA15I across the boundary, i.e. from — e to + e. From this integration we obtain 
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Inserting the solutions given in equations <A19> in the boundary conditions results in four linear equations for A\, B\, A2 and 62- 
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Here we have used 
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and 
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in order to preserve the symmetry inherent in the boundary conditions. The set of linear equations can be conveniently written in matrix form, 
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where we used the definitions of q\ and q 2 , divided all elements in the fourth row by (pj + p 2 ) and set 
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The determinant of the matrix must vanish and this condition is equivalent to the required equation for n. We simplify the calculation of the 
determinant by dividing the elements of the first and second column by n\ as well as the elements of the third and fourth column by n' 2 . We 
then subtract the first column off the second and the third column off the fourth. We then multiply the elements of the first column by n' x /n' 2 
and add this column to the third column. Finally, we multiply the elements in the third column by n' 2 and, after again using the definitions of 
q\ and q 2 , we obtain the reduced determinant 

k fn'j + « 2 ) q_2 — k 



?1 



a in. 



2k 2 (aiVin'j -a 2 v 2 n 2 ) 



&\n\ +C(l- 



I 1 (i - ^r) t\—\. («, —11, 1 ~ t«j«; ; 

The equation resulting from expanding the above determinant is complicated. We do not reproduce it here, but note that it correctly reduces 
to equation (1 13) in chapter X of C61 for the case of U\ — U 2 — 0. 

In this paper we are specifically interested in cases where the mass density of the heavy fluid is much greater than that of the light fluid, 
i.e. p 2 S> pi. In this case, ai ~ and a 2 ~ 1 and expanding the determinant in equation IA31i yields 
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R + Ciri, 
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-a 2 ri 2 2 ) 



-a 2 n 2 



a 2 n 2 — C (l - 
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-0. 



(A31) 



q 2 ' + 2k 2 ql-4k i q 2 



0, 



(A32) 



where we substituted n' 2 = (q 2 , — k 2 ) v 2 . To arrive at equation JA32I we have already divided out the trivial solutions n' 2 = 0, q\ — k = and 
k 2 — q 2 = 0, all of which only describe oscillations of the boundary with constant amplitudes. 

Several interesting conclusions can be drawn from our results. The properties of the light fluid are not important in determining the 
behaviour of instabilities at the boundary between the two fluids, since equation <A32> does not depend on any quantities with a subscript 
'1'. Also, the velocities U\ and U 2 do not explicitly appear in equation <A31> . Let A'] be an arbitrary solution of equation <A31> for n' x and 
N 2 the corresponding solution for n' 2 . Then 



n = N\ - ik x ll\ 
n = N 2 — ik x U 2 



(A33) 



give possible growth rates for instabilities. Since A'i and N 2 must be independent of U\ and U 2 , these two velocities can only give rise to 
oscillations of the boundary with fixed amplitudes. In other words, in the linear regime K-H instabilities cannot grow on the boundary between 
two viscous fluids. This result is due to the last term in equation <A15> . which requires the continuity of dw/dz for non-zero viscosities fj. In 
the inviscid case K-H instabilities can grow, but only for wavenumbers exceeding (C61) 



(A34) 



aia 2 (£/i-£/2r 

For our case of pi <C p 2 we have fc m ; n — > °° and so the K-H instabilities are suppressed on all length scales. 

The 4 th -order polynomial in q 2 of equation <A32t has four analytical, in general complex, solutions. The expressions for these solutions 
are very complicated and, unfortunately, cannot be turned into simple prescriptions for the existence of solutions for n with positive real 
parts. However, in this paper we present plots of these solutions for the specific conditions investigated. 

If both fluids are inviscid, we note that v 2 q 2 = and we recover the solution presented in C61, i.e. 



Igk 1 



gp2/' 



(A35) 
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In this form it is straightforward to see that surface tension leads to a. maximum wavenumber £max — VsPl/Ts, so that instabilities with 
k > k max cannot grow (C61). The same fc max also applies in the viscous case. 

Finally, in the case of inviscid fluids without surface tension at the boundary we recover the well-known result 




